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Abstract. The unitary group U(N) acts by conjugations on the space Jf?(N) 
of A x A Hermitian matrices, and every orbit of this action carries a unique 
invariant probability measure called an orbital measure. Consider the projection 
of the space Jf?(N) onto the real line assigning to an Hermitian matrix its (1, 1)- 
entry. Under this projection, the density of the pushforward of a generic orbital 
measure is a spline function with A knots. This fact was pointed out by Andrei 
Okounkov in 1996, and the goal of the paper is to propose a multidimensional 
generalization. Namely, it turns out that if instead of the (1, l)-entry we cut out 
the upper left matrix corner of arbitrary size A x A, where A = 2, . . . , A — 1, 
then the pushforward of a generic orbital measure is still computable: its density is 
given by a A x A determinant composed from one-dimensional splines. The result 
can also be reformulated in terms of projections of the Gelfand-Tsetlin polytopes. 



Contents 



1. Introduction 

2. Preliminaries 

3. Projections of orbital measures 

4. Acknowledgement 
References 



1. Introduction 

Orbital measures. Let Jff(N) be the space oi N x N Hermitian matrices. For 
K = 1,...,N — 1, we denote by p 1 ^ : Jif(N) — > J(?(K) the linear projection 
consisting in deleting from the matrix H e Jif(N) its last N — K rows and columns. 
We call Px(H), the image of H under this projection, the K x K corner of H. 

The unitary group U(N) acts on Jf(N) by conjugations, and because U(N) is 
compact, each orbit of this action carries a unique invariant probability measure, 
which we call the orbital measure. Given an orbital measure \x on Jif(N), denote 
by Px([i) its pushforward under projection p^. Our goal is to describe p^jj). 

The orbits in Jrff(N) (and hence the orbital measures) can be indexed by A^-tuples 
of weakly increasing real numbers X = (x\ < ■ ■ ■ < x^), the matrix eigenvalues. 

l 
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Let &(N) C R N denote the set of all such X's. Given X G 3£{N), we write O x 
and /xx for the corresponding orbit and orbital measure, respectively. 

Since p^i^x) is a £/(i^)-invariant probability measure on J$?(K), it can be uniquely 
decomposed into a continual convex combination of orbital measures, governed by 
a probability measure vx,k on the parameter space 2£{K). That is, vx,k is charac- 
terized by the property that, for an arbitrary Borel subset S C JT(iV), 

(PkM)(S) = [ Hr{S)v x , K {dY). 

JYeX(K) 

The measure Ux,k can be called the radial part of measure Pxinx)- 

Main result. Denote by JT°(iV) the interior of <5T(iV); that is, i*T°(iV) consists of 
iV-tuples of strictly increasing real numbers. If X G 3£ Q (N), then is absolutely 
continuous with respect to Lebesgue measure on 3£(K) C M. K , and the main result, 
Theorem I3.3[ gives an explicit formula for the density of vx,k- 

In the case K = 1 the target space of the projection is the real line, and the density 
in question coincides with a B-spline, a certain piecewise polynomial function on R 
(this fact was observed by Andrei Okounkov). In the general case, it turns out that 
the density of vx,k is expressed through a K x K determinant composed from some 
B-splines. 

As the reader will see, the proof of Theorem l3.3l is straightforward and elementary. 
The main reason why I believe the result may be of interest is the very appearance 
of splines, which are objects of classical and numerical analysis, in a problem of 
representation-theoretic origin. 

Gelfand— Tsetlin polytopes. Before explaining a connection with representation 
theory I want to give a different interpretation of the measure vx,k- 

For X G 3£(N) and Y G 3£{N — 1), write Y -< X or X y Y if the coordinates 
of X and Y interlace, that is 

xi < yi < x 2 < ■ ■ ■ < x N ^ < y N ^x < x N . 

Given X G 3£{N), the corresponding Gelfand- Tsetlin polytope Px is the compact 
convex subset in the vector space 

R N-1 x R N-2 x ... xR = E iV(iV-l)/2 ) 

formed by triangular arrays subject to the interlacement constraints: 

P x : = {(F^-D, . . . , FW) e R N( - N -^ 2 : X y Y( N -V yd)}. 

Consider the map assigning to a matrix H G Ox the array formed by the collections 
of eigenvalues of its corners p^_ 1 (H),p^_ 2 (H), . . . ,p^(H). It is well known (see 
Corollary 13.21 below) that this map projects the orbit Ox onto the polytope P x 
and takes [l x to the uniform measure on Px (that is, the normalized Lebesgue 
measure). Next, given K = 1, . . ., N — 1, consider the natural projection Px — > 
3£{K) extracting from the array (Y^-V, yd)) its Kth component Y^ K \ The 
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measure vx,k is nothing else than the pushforward of the uniform measure under 
the latter projection. 

Discrete version of the problem: relative dimension in Gelfand— Tsetlin 
graph. Let G1V := ^(N) fl 7L N be the set of weakly increasing iV-tuples of 
integers. The elements of GT^ are in bijection with the irreducible representations 
of the group U(N): with X = (xx, . . . , Xjy) G GT N we associate the irreducible 
representation T x with signature (=highest weight) X := (xn, . . . , x\). Here we 
pass from X to X, because the coordinates of signatures are usually written in the 
descending order, see Weyl [13]. 

Let X G GTn and consider the finite set Pf := Px H Z N<yN ~ 1 ^ 2 consisting of 
integral points in the polytope Px- Let us replace the uniform measure on Px by 
the uniform measure on Pf (that is, the normalized counting measure). Next, given 
K — 1, . . . , N — 1, we consider again the same projection Px — > S£{K) as before 
and denote by v x K the pushforward of the uniform measure on P\ . Evidently, v\ K 
is a probability measure with finite support. 

Elements of P\ are the Gelfand-Tsetlin schemes (also called Gelfand-Tsetlin 
patterns) with top row X; they parameterize the elements of Gelfand-Tsetlin basis 
in Tx- By the very definition of v\ K , for Y G GT^, the quantity v\k^X) (the mass 
assigning by v\ K to Y) equals the fraction of the schemes with the Kth row equal 
to Y. This quantity is the same as the relative dimension of the isotypic component 
of Ty in the restriction of T x to the subgroup U(K) C U(N). 

The Gelfand-Tsetlin graph has the vertex set GTi U GT 2 U . . . and the edges 
formed by couples Y -< X. In the terminology of Borodin-Olshanski [3], v x K (Y) is 
the relative dimension of the vertex Y G GT^ with respect to the vertex X G GUV. 
In [5], we derived a determinantal formula for the relative dimension (see also Petrov 
[TTJ for a different proof). That formula can be viewed as a discrete version of the 
formula of Theorem 13.31 

I first guessed the formula of Theorem 13.31 by degenerating the "discrete" formula 
of [3]. However, this is not an optimal way of derivation, because the discrete case 
is much more difficult than the continuous one. I am grateful to Alexei Borodin for 
the suggestion to study the degeneration of the "discrete" formula. Note that from 
the comparison of the measures vn,k an d Vx k ^ ^ s seen ^ na ^ the former should be 
related to the latter by a scaling limit transition. 



2. Preliminaries 

The fundamental spline with n > 2 knots yx < ■ ■ ■ < y n can be characterized 
as the only function a h- > M(a; y\ . . . , y n ) on R of class C n ~ 3 , vanishing outside the 
interval (yx, y n ), equal to a polynomial of degree < n — 2 on each interval (yi, yi+%), 
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and normalized by the condition 

X) 

M(a;y 1 , . . .,y n )da = 1. 



Here is an explicit expression: 



M(a; y u . . . , y n ) := (n - 1) £ ^ ^ ■ (1) 



*«>. n («•-!/■ 

In particular, for n = 2 



M(a ;2/l ,y 2 ) = ^i^. 

2/2 -yi 

Remark 2.1. The above definition is taken from Curry-Schoenberg In the 
subsequent publications, Schoenberg changed the term to B-spline. The latter term 
became commonly used. However, in the modern literature, it more often refers to 
the function 

B(a; yi,...,y n ):= (y n ~ Vl) Y] Jr , ^ r , (2) 

which differs from M(a; yi, . . . , y n ) by the numerical factor (y n — y\)j (n — 1); see, 
e.g., de Boor [2] or Phillips [12] . The normalization in (j2J) has its own advantages, 
but we will not use it. Note also that M(a; j/j, . . . , y n ) is a special case of Peano 
kernel, see Davis [5], Faraut [7]. 

We need two well-known formulas relating M(a; yi, . . . , y ra ) to divided differences 
(see, e.g., @], 0). 

Recall that the divided difference of a function f(x) on points yi, . . . , y n is defined 
recursively by 

/ 2/1,3/2 = , / 2/i> 2/2,2/3 = , 

2/2-2/1 2/3 - 2/1 

and so on; the final step is 

,r 1 f[V2, ■ ■ ■ ,2/rJ - f[Vl, ■ ■ ■ ,2/n-l] 

/2/l,---,2/n = • 3 

Vn - 2/1 

Next, set 

x s , x > 
0, x < 0. 

In this notation, the first formula in question is 

M(a;y 1 ,...,y n ) = (n-l)f[y 1 ,...,y n ], where f(x) := (x - a)™~ 2 , (4) 



x + 



and the second formula is 

/[</!,..., y n ] = 7^1)? / M(a;y 1 ,...,y n )f^- 1 \a)da. (5) 

In (jSJ), / is assumed being a function on R with piecewise continuous derivative 
of order n — 1. In particular, 03) is applicable to f(x) — (x — £)+ , which is used in 
the lemma below. 

To shorten the notation, let us abbreviate Y := (y± < ■ ■ ■ < y n ). 

Lemma 2.2. Fix n = 2, 3, . . . and an n- tuple Y = (y 1 < ■ ■■ <y n ) e S£ Q {N). For 
an arbitrary b el set 

f b {x):={x-b) n + -\ xeR. 

One has 

! M(a;Y)da=l-f c [Y], c G R, (6) 

J —oo 

J C M(a;Y)da = f b [Y]-f c [Y], b < c, (7) 

/+oo 
M(a;Y)da = f b [Y], beR. (8) 

Proof. To check ((?]), we apply (jSD to /(x) = — / c (x), which is justified, see the 
comment just after (jSJ). Then in the left-hand side of (jSJ) we get — / c [^]- Next, 
observe that the (n — l)th derivative of equals (n — l)\l x >b, so that 

/(«-i)(a) = ( n - l)!(l a > 6 - l a > c ) = (n - 1)!1 6 < 0<C . 

Therefore, in the right-hand side we get J fe ° M(a; y)rfa, which proves (j7j). 

Now (|S]) follows from ^7]) by setting c = +oo, and (jSJ) follows from (jSJ), because 
the total integral of the M(a; Y) equals 1. □ 

3. Projections of orbital measures 
We keep to the notation of Sections [TJ and [2] 

Given X G 3£(N), the pushforward of the orbital measure \ix under the map 

Ox 3 H \-¥ the spectrum of p%-\(H) 

can be viewed as a probability measure on JT(iV— 1) depending on X as a parameter; 
let us denote it by k%_ x {X, ■ ) or A^_ X (X, dY). We regard A(X,dY) as a Markov 
kernel. 

By classical Rayleigh's theorem, the eigenvalues of a matrix H G J%~(N) and its 
corner interlace. Therefore, the measure Aj^_ 1 (X, ■ ) is concentrated on 

the subset 

{Y g 3£ (N - 1) : Y -< X} c R N -\ (9) 
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Proposition 3.1. Assume X = (xi, . . . ,xn) £ =3T (iV). Then the measure A ^^(X, 
is absolutely continuous with respect to Lebesgue measure on the set and the 
density ofA^_ ± (X, ■), denoted by Aj^_ 1 (X, Y), is given by 



A^_ 1 (X,F) = (iV-l)!^|||w, (10) 



where we use the notation 

V(X) = J]( 

j>i 

and the symbol 1y-<x equals 1 or depending on whether Y -< X or not. 

Proof. To the best of my knowledge, a published proof first appeared in Baryshnikov 
PQ Proposition 4.2]. However, the argument given in [1] was known earlier: it is 
hidden in the first computation of the spherical functions of the groups SL(N, C) 
due to Gelfand and Naimark, see [8j §9]. Note also that a more general result can 
be found in Neretin [9]. 

Here is a different proof. Consider the Laplace transform of the orbital measure 

£x(Z) := J e T ^ ZH ^ x (dH), (11) 

where Z is a complex N x N matrix. The integral in the right-hand side is often 
called the Harish-Chandra-Itzykson-Zuber integral. Its value is given by a well- 
known formula (see, e.g., Olshanski-Vershik [TOj Corollary 5.2]): 

fl x (Z) = c N J,J ~ 1 12 

{[{Zj - z i )(x j - Xi) 
j>i 

where z\ , . . . , z x are the eigenvalues of Z and 

c N = (N- 1)\(N- 2)!... 0! 

(note that the right-hand side of ( jl2j) does not depend on the enumeration of the 
eigenvalues of Z). 

The claim of the proposition is equivalent to the following equality: Assume that 
the entries in the last row and column of Z equal 0, so that Z has the form 



Z 




(13) 

where Z is a complex matrix of size (N — 1) x (N — 1); then 

MZ) = J V(Y)MZ)dY. (14) 

To prove (Tl4"|) . consider the matrix T := [e ZlXj ] in the right-hand side of ffl2l) . 
Since Z has the form ( fl3l) . at least one of the eigenvalues Z\, . . . , zn equals 0. It is 
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convenient to slightly change the enumeration and denote the eigenvalues as zq = 
0, Zi, . . . , Zn-x- In accordance to this we will assume that the row number i of T 
ranges over {0, . . . ,N — 1} while the column index j ranges over {1, . . . , N}. Since 
Zq = 0, the 0th row of T is (1, . . . , 1). Let us subtract the (N — l)th column from 
the iVth one, then subtract the (N — 2)th column from the (N — l)th one, etc. This 
gives det T = det T, where T stands for the matrix of order N — 1 with the entries 

%. = e **i+i _ e **i = Zi e*Vd yj , i,j = l,...,N-l. 



It follows 



so that 



detT = Zl ...z N [ dY det[e z ^}?- = \, 

JY<X 

zx...z N f Y ^ x dYdet[e^}»- = \ 



fix(Z) = c N 



n - ^) ■ v(x) 

N-l>j>i>0 



Next, because z = 0, the product over j > i in the denominator equals 

Z\ . . . Zn J~J (zj — zi), 

N-l>j>i>l 

so that the product Zi . . . z N is cancelled out. Taking into account the fact that 
is given by the determinantal formula similar to (fl2|) and using the obvious 
relation cn = (N — l)!c/v-i we finally get the desired equality (JHJ). 

□ 

From Proposition 13.11 it is easy to deduce the following corollary (see also (TJ 
Proposition 4.7] and Defosseux [6]). 

Corollary 3.2. Fix X e JT°(iV) and let H range over Ox- The map assigning to H 
the collection of the eigenvalues of the corners p^(H) , where K = N—l, N—2, . . . , 1, 
projects Ox onto the Gelfand-Tsetlin polytope Px and takes the measure fix to the 
Lebesgue measure multiplied by the constant 

(N — 1)1 (N — 2)! ... 0! 

V{X) ' 

In particular, the volume of Px in the natural coordinates is equal to the inverse of 
the above quantity. 

Recall that vx,k stands for the radial part of the K x K corner of the random ma- 
trix H G Ox, driven by the orbital measure \ix (see Section [TJ), and Mia; y x , . . . , y n ) 
denotes the fundamental spline with n knots yi, . . . ,y n (see ([I]) and (j3J)). 
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Theorem 3.3. Fix X = (xi, ...,x N ) G Jf°(A). For any K = 1, . . . , N - 1, the 

measure vx,k on $f(K) is absolutely continuous with respect to Lebesgue measure 
and has the density 

V(A) det [M(a j ;x i ,...,x N - K+i )]« =1 

M [a,i, . . . , a K , Xx,...,x N ) .— c N)K — — — r , [io) 

11 \ x j x i) 

j-i>N-K+l 

where 

j-r f N — K + i 
cn,k =11 

8=1 



I 



Note that for K = 1 the right-hand side reduces to the fundamental spline with 
knots xi, . . . ,xn- Thus, in the case K — 1 the theorem says that the density of 
the measure u N> i on R coincides with the spline M(a; x\, . . . ,xjy). This simple but 
important claim is due to Andrei Okounkov, see [TO] Proposition 8.2]. 

Proof. We argue by induction on K, starting with K = N — 1 and ending at K = 1. 

Step 1. Examine the case K = N — 1, which is the base of induction. We 
have ux,N-i{dA) = A^_ 1 (X, dA). By proposition 13. 1\ the measure A]^_ 1 (X, •) on 
JT(A— 1) is absolutely continuous with respect to Lebesgue measure and has density 
A^_ t (X, A) given by ({TO]) . Thus, we have to check that A^_ 1 (A, A) coincides with 
the quantity M(ai, . . . , ajy-i; %i, ■ ■ ■ , Xn) given by the right-hand side of (TT5|) . where 
we have to take K = N — 1. That is, the desired equality has the form 

,„ ,„V(A), V{A)Aet[M(a j ;x i ,x i+1 )\l-} 1 
(N - 1)1^ W = n . 

Since cjy^-i = (A — 1)!, the desired equality reduces to 
det [M(aj] x h x i+1 )]^~J l 



(x 2 - X 1 )(X 3 ~X 2 ) ... (X N - Xn-x) ' 

Observe that the (i,j)-entry in the determinant is the quantity 

M{a j ;x i ,xu. 1 ) = , 

X%-\-\ Xi 

which vanishes unless aj G [x^Xj+i]. Since a\ < ■■■ < on-i, the determinant 
vanishes unless A -< X. Furthermore, if A -< X, then the matrix under the sign 
of determinant is diagonal, so the determinant equals the product of the diagonal 
entries, which equals 

1 



(x 2 - x 1 )(x 3 -x 2 ) ... (x N - X N . 

as required. 



Step 2. Given K = 1, . . . , N — 1, we consider the superposition of Markov kernels 

\N ._ aV aTV-1 

ly K ■— IV N-1 1V N~2 ■ ■ ■ Ly K ■ 

In more detail, the result is a Markov kernel on JT(iV) x 3£{K) given by 
A%(X, dA) = [a^X, dY^A^Y^^Y^) . . . A K K + \Y^ K+1 \ dA), 



where the integral is taken over variables y^" 1 ), . . . , Y^ K+1 \ Obviously, A%(X, dA) = 
vx,i<{dA), which entails the recurrence relation 

vx,k-x = vx,kAI_ x , K = N-l,N-2,...,2, (16) 
where, by definition, v XK A^_ x is the measure on S£{K — 1) given by 



{v x , K A K K _ x ){dB) = / u x>K (dA)A^_ 1 (A,dB). (17) 

Step 3. Assume now that the claim of the theorem holds for some K > 2 and 
deduce from this that it also holds for K — 1. To do this we employ ( II 6p and ( II 7J . 

First of all, (1T61) and (fl71) imply that Px,k-i is absolutely continuous with respect 
to Lebesgue measure on S£{K — 1) and has the density 

(y X)K kK_ x ){B) = [ v XjK {dA)A^_ 1 {A, B), Be ^(K - 1). (18) 

Let us compute the integral explicitly. By the induction assumption, v Xt K is 
absolutely continuous and has density ( IT5l) . Therefore, integral ( !T8l) can be written 
in the form 

/ M(ai, . . . , a K \ Xi, . . . , xtv)A^_ 1 (A, B)da\ . . . dax- 

Write B = {pi, . . . ,&&-). Substituting the explicit expression for A^_ 1 (A,_B) given 
by Proposition 13.11 we rewrite this as 

(K - 1)\V(B) J M( ai ,...,a K ^...,x N ) dai ^ (ig) 

where the integration domain is 

- oo < ai < bx, h < a i+1 < b i+ i, 6 fc -i < cl k < +oo. (20) 

Next, plug in into ( TT9|) the explicit expression for M(a%, . . . , a K ; x±, . . . , x^r) given 
by ({15]) . Then the factor V(A) is cancelled out and we get 

c N>K (K-l)\V(B) f K 
— — - — r / det[M(a j ;x i ,...,x N -K+i)\ id=1 dai...da K (21) 

11 \ x j x i) J A 

(j,i):j-i>N-K+l 

with the same integration domain (!20l . 
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Put aside the pre-integral factor in ( 12TT) and examine the integral itself. It can be 
written as a K x K determinant, 

det[F(z,j)]£ =1 , 

where 



F(i,j):= / M^a-Y^da 

and 

with the understanding that b = — oo and bx = +00. 
We are going to prove that 



det[F(i, = (N-K+ l)^" 1 H {x N . K+l+l - 



K-l 

...... 

8=1 



x det[M(6 i ; x i} . . . , x N „ K+i+1 )]«r} v (22) 
This will justify the induction step, because 

(N-K+l) K - 1 



cn,k — CN,K-1 



{K-l)\ 
and 

Efi=l {xN-K+i+l - Xi) _ 1 



n - Xi) n i x j - x i) ' 

j-i>N-K+l (j,i):j-i>N-K+2 

Step 4. It remains to prove (|22|) . We evaluate the quantities F(i,j) using Lemma 



I2.2[ where we substitute n = N — K + 1 and Y — Yi. Then we get that the matrix 
entries F(i,j) are given by the following formulas: 

• The entries of the first column have the form 

F(z,l) = l-/ 6l [^] by®. 

• The entries of the jth column, 2 < j < K — 1, have the form 

F(i,j) = f h ._M-f bj \Y i ] by®. 

• The entries of the last column have the form 

F{i,K) = f bK [Y i \ by®. 
We have det F = det G, where the K x K matrix G is defined by 

G(i,j) :=F(i,j) + --- + F(i,K). 
The entries of the matrix G are 

G(i,l) = l, G(i,j) = f bj _ M, 2<j<K. 
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Next, we get det G = det H with the (K — 1) x (K — 1) matrix H defined by 
H(i,j) := F(i + 1, j + 1) - F(i,j), \<i,j<K-l. 
Observe now that 

H{iij) = f h .[Y i+1 \-f b ,[Yl 

which can be rewritten as 

H(i,j) = (x N -K+i+i - Xi)— - (23) 

XN-K+i+l " x i 

= (x N -K+i+i ~ Xi)f bj [xi, x N -K+i+i] by §3§ (24) 

= Tt \? , i ( x N-K+i+i - Xi)M(bj; Xi,..., x N _ K+i+ i) by (0}. (25) 

iv — K + 1 

This shows that the determinant det H = det[.H"(i, i)]^-=i equals the right-hand side 
of (T2"2l . Since det H = det F, this completes the proof. 

□ 
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